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EIG (Compute Eigenvalues and Eigenvectors) 


DECUS Program Library Write-up 


DECUS No. 8-317 


ABSTRACT 

The extraction of eigenvalues and eigenvectors is carried out by means of a group of pro¬ 
grams, the intermediate stages of the computations being stored on the disk. The eigen¬ 
values, and the associated eigenvectors, are extracted one by one, starting with the 
largest eigenvalue, and the process can be terminated as soon as sufficient eigenvalues 
have been extracted from the basic data matrix. The programs were designed to enable 
principal component and canonical variate analyses to be carried out. 

TAPES REQUIRED 


1 • Form of program tapes - All program tapes are written in the PDP-8 FORTRAN-D 
language, and are in the source language. There are four programs in the set, as follows: 

(a) SMOD - program to store data matrix on disk. 

(b) EIG1 - program to compute the highest latent root and vector. 

(c) EIG2 - program to extract and print the highest latent root and vector. 

(d) AITK - program to square the current data matrix, so as to increase the speed of con¬ 
vergence in the process of extracting the latent root. 

2* Form of data tape - The data tape should consist of the matrix of correlation coefficients 
or of corrected sums of squares and products. The program requires the complete matrix, 
including the principal diagonal. 

OPERATING INSTRUCTIONS 

Operation of this suite of programs is simplified if all of the programs are compiled onto 
the disk before the data are read, i. e. 

.FORT 

*OUT-S:SMOD 

* 

*IN-R: 

* 

*. FORT 

*OUT-S:EIGl 

* 

*IN-R: 

* ''p* 











*. FORT 

*OUT-S:EIG2 

* 

*IN-R: 

* T 

*. FORT 

*OUT-S:AITK 

* 

* IN-R: 

* T 


The data tape should now be placed in the high-speed reader, and read onto the disk as 
follows: 


. FOSL 

*IN-S:SMOD 


*OPT-S: 

*OUT-S:DATA 

*IN- 

* 

* READY 
T 


The program will pause for the entry of the number of variables, i. e. the dimension of the 
matrix, and this should be entered on the teletype and terminated by "Return. " The program 
will type on completion. 

The first latent root and vector can now be computed, as follows: 


.FOSL 

*IN-S:EIG1 


*OPT-S: 

*OUT- 

*IN-S:DATA 

* 

*READY 

T 


The program will pause for the entry of the dimension of the matrix, and this should be 
entered on the teletype and terminated by "Return. " The high-speed punch should also 
be switched on as soon as this has been done. The program will output a vector on the 
high-speed punch and then type on completion. The speed of convergence of the com 
putation depends partly on the dimension of the matrix and partly on the closeness of the 
latent roots. If the calculation does not end after five minutes or so, it is advisable to re 
start the computation after squaring the matrix (see operating instructions for the AITK 
program). 
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The latent root and vector can be printed, and extracted from the data matrix, as follows: 

. FOSL 

*IN-S:EIG2 

* 

*OPT-S: 

*OUT-S:DATA 

*IN-S:DATA 

Output from EIG1 in high-speed 
*READY reader 

T 

The program will pause for the entry of the dimension of the matrix, and this should be 
entered on the teletype and terminated by "Return." The program will print the eigenvalue 
and the eigenvector, and then store the residual data matrix on the disk, typing 1 on com¬ 
pletion. If further latent roots and vectors are required, return to program EIG1 above. 

Convergence of the iterative procedure used in calculating the latent root and vector may, 
in certain cases, be slow. In these cases, the data matrix can be modified to improve 
the convergence as follows: 

.FOSL 

*IN-S:AITK 

* 

*OPT-S: 

*OUT-S:DATA 

*IN-S:DATA 

* r 

*READY 

r 

The program will pause for the entry of the dimension of the matrix, and this should be 
entered on the teletype, and terminated by "Return." This procedure can be repeated 
three or four times. In general, the speed of convergence will be increased by a factor 
of 2 for every time the procedure is used, but it is not advisable to apply the procedure 
too often because of possible loss of accuracy. 

In difficult cases, convergence may not be reached even by applying the Aitken process. 

In such cases, the condition for convergence may be relaxed by altering the FORTRAN 
statement following the labelled statement 9 to: 

IF (TEST-0.001) 6,6,7 


OUTPUT 

The only printed output is obtained from the EIG2 program, ond consists of the eigenvolue 
and the associated eigenvector. 

STORAGE 

Normal for FORTRAN-D 

For the 4K version, the dimension of the matrix is limited to 12. 








METHOD 


See "A Short Course In Multivariate Analysis/ 1 M. G. Kendall, Griffin, 1956. 


SMOD 


L 

C PROGRAM TO STORE MATRIX ON DISK 

DEFINE DISK 
DIMENSION R(144) 

ACCEPT 100 , N 
N=N*N 
DO 1 1=1, N 
READ 2,101,R(I) 

1 CONTINUE 
DO 2 1=1, N 
WRITE 3,lj01,R(l) 

2 CONTINUE 

100 FORMAT (I) 

101 FORMAT (/,E) 

END 









EIG 1 


L 

C 


100 

101 



1 

102 


2 


3 

7 


4 


5 


9 

6 
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PROGRAM TO COMPUTE LATENT ROOT AND VECTOR 
DEFINE DISK 

DIMENSION R(144),A(12),B(12) 

TYPE 100 

FORMAT (/, "ENTER NO OF VARIABLES",/,/) 

ACCEPT 101, N 
FORMAT (I) 

NN=N*N 
DO 1 1=1 ,NN 
READ 3,102, R(l) 

CONTINUE 
FORMAT (/, E) 

DO 2 1=1, N 
A(l)=0.0 
DO 2 J=1, N 
K=I+N*(J-1) 

A(I)=A(I)+R(K) 

CONTINUE 

V=A(1) 

DO 3 1=1, N 

A(I)=A(I)/V 

CONTINUE 

TEST=0.0 

K=1 

DO 4 1=1, N 
B (l)=0.0 
CONTINUE 
DO 5 1=1, N 
DO 5 J=1, N 
B(J)=B(J) +R(K)*A(I) 

K=K+1 

CONTINUE 

V=B(1) 

DO 9 1=1, N 
B(I)=B(I)/V 

TEST=TEST+ABSF(A(I)-B(I) ) 

A(I)=B(I) 

CONTINUE 
IF (TEST-0.0001) 6,6,7 
DO 8 1=1, N 
WRITE 2,102, A(l) 

CONTINUE 

END 
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EIG 2 


L 

C 


100 

1 

101 


2 


3 


4 


5 


6 


7 


PROGRAM TO PRINT LATENT ROOT AND VECTOR 
DEFINE DISK 
DIMENSION R(144),A(12) 

ACCEPT 100, N 
FORMAT (I) 

DO 1 1=1, N 
READ 2,101 ,A(I) 

CONTINUE 
FORMAT (/, E) 

NN=N*N 
DO 2 1=1, NN 
READ 3, 101, R(l) 

CONTINUE 
ROOT=0.0 
DO 3 1=1, N 
ROOT=ROOT+A(l)*A(l) 

CONTINUE 

ROOT=l .0 /SQTF(ROOT) 

EIG=0.0 
DO 4 1=1, N 
EIG=EIG+A(I)*R(I) 

CONTINUE 
TYPE 101, EIG 
DO 5 1=1, N 
A(l)=ROOT*A(l) 

TYPE 101, A(l) 

CONTINUE 
DO 6 1=1, N 
DO 6 J=1, N 
K=I+N*(J-1) 

R(K)=R(K)-A(I)*A(J)*EIG 
CONTINUE 
DO 7 1=1, NN 
WRITE 3, 101, R(l) 

CONTINUE 

END 


* 
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AITK 


PROGRAM TO SQUARE MATRIX 
DEFINE DISK 
DIMENSION R(144) 

ACCEPT 100, N 
FORMAT (I) 

NN=N*N 
DO 1 1=1, NN 
READ 3,Ufl,R(|) 

CONTINUE 

IR=0 

IK=-N 

DO 10 K=1,N 

IK=IK+N 

DO 10 J=l, N 

IR=IR+1 

JNJ-N 

IB=IK 

RNEW=,0.j0 
DO 20 1=1, N 
JNJI+N 
IB=IB+1 

RNEW=RNEW+R(J|)*R(|B) 
CONTINUE 
WRITE 3,101, R NEW 
CONTINUE 
FORMAT (E) 

END 
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